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7 OEM TJX AS TOE USB IB BOUHDAEY- LATHE CALCTTLATIOHS 

OH LOW-DEAO WIH&S 
By 1. IT. Jacobs and A. I. von Doenhoff 



IHTEODUOTIOH 



This report has - been prepared in response to frequent 
requests received from aircraft companies for information 
concerning methods of computing on low-drag wings the 

(l) Transition point 

(3) Telocity distribution in laminar and turbulent 
boundary layers 

(3) Thickness of the boundary layer, both laminar 
and turbulent 



Location of Transition Point 

The best estimate of the transition point that we 
can make at the present time is based on the data obtained 
from tests of a glove on the B-*18 airplane in flight* 
Under one set of conditions values of Eg in the laminar 
boundary layer were observed between 8000 and 9500, where 
Eg b V6/v and 7 is the velocity outside the boundary 
layer at the point in question, 8 the distance normal 
to ttie surface from the surface to the point in the bound- 
ary layer where the velocity head is one-half its local 
value outside the boundary layer (u « 0.707 V; see fig* 
l) , v the kinematic viscosity. On the new low-drag 
sections the laminar separation point is so close behind 
the point of minimum pressure that transition is assumed 
to occur at the minimum pressure point if E. is less 
than the previously mentioned value at the minimum pres- 
sure point* 

"rfe have not made any transition measurements on 
smooth streamline bodies under completely turbulence-free 
conditions. It Beems reasonable to suppose, however, that 
the oondition for transition is not very different from 
that for two-dimensional bodies. 
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Velocity Distribution in Laminar and 



Turbulent Boundary Layers 



Boundary-layer velocity distribution measurements 
have shown that in a favorable pressure gradient the 
shape of the laminar distribution is very closely approxi- 
mated by the Blasius distribution for a flat plate. An 
approximate relation for this distribution! due to Pohl- 
hausen, is given by the formula 



The outer limit of the boundary layer corresponds to 



so soon In a region of adverse pressure gradient that 
this region is of little importance. The velocity distri- 
butions in this region can, however, be calculated by the 
methods described in reference 1, although this is usually 
a very lengthy process. 

In a turbulent boundary layer, the 1/7-power law 1b 
a fair approximation to the shape of the velocity profile. 



favorable. In an unfavorable gradient there is a gradual 
departure from this shape, the departure increasing with 
the distance behind the point of minimum pressure. The 
problem of finding the shape of the turbulent boundary- 
layer velocity distribution as a function of pressure 
gradient has not yet been solved. 



Thickness of Laminar and Turbulent Boundary Layers 

The thickness of the laminar layer in a region of 
favorable pressure gradients is oomputed by assuming that 
the shape of the velocity distribution is the same as the 
Blasius distribution and then integrating the Zarman 
momentum relation. The result may be given in the fol- 
lowing form for two-dimensional flow 





In most cases laminar separation occurs 




when the pressure gradient is 
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H c » — ^— wine Heynolds nujaber 

s distance along surface from leading edge 
c chord 

7 Q reference velocity, usually token bb the free 
stream velocity 

7- velocity outside boundary layer at any value 
of s / c 

7j velocity outsiflo "boundary layer at point for 
which "boundary layer is "being conputed 

The corresponding relation for bodies of revolution 
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where Hj, = 



L length of "body 

r distance of surface from axis of revolution 
at any value of s/^i 

r x distance of surface from axis of revolution 

at point for which "boundary layer is being 
computed 
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The most reliable estimates of the momentum thickness 
of the t\irbulent boundary layer in t\vo-dimensional flow 
may be found from the relations given in reference 2. 
These relations, which are valid in those cases where 
there is no danger of turbulent separation,' are as follows S 
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The quantity G is the momentum thickness of the 

boundary layer 6 - -=- / u(V-u)dy. lor the Blasius 

•Jo 

distribution G = 0.289 8. The initial value of 9 for 
the turbulent layer is taken the same as that for the 
laminar boundary layer at the transition point. The cor- 
responding value of £ for the turbulent boundary layer 
is then found from equation (l). The variation of £, 
and hence by eqiiation (i) of 6, along the surface is 
then found from equation (2), which is solved by a step- 
by-step process. In more detail, if t n is the known 

value of £ at any point s/c along the surface, 
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where t is the value of t at (s/c) + A (s/c). The 

process is then repeated, using £ r+1 . as the- new value 
of £ n . The value of Ms/c) is arbitrary within limits. 
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If — — is not varying rapidly with s/c, A(s/c) 

* (f ) 



on may be taken larger than when rapid variations in 



(I) 



occur. In most computations the value of Ms/c) ma y 
be taken as 0.01. 

The quantity £ is related to a skin friction coef- 
ficient 

t =~ 

where T shearing stress at surface 
p dens i ty 

The profile drag of the airfoil can be obtained from 
the values of 0 at the trailing edge on the upper and 
lower surfaces from the following equation 



_ gCgg + 9l) / Vt\ 5,s 
c o 



where 9^ value of 8 at the trailing edge, upper sur- 
face 

8 L value of 9 at the trailing edge, lower sur- 
face 

Vj. velocity outside the boundary layer at the 
trailing edge 

For bodies of revolution the relations given in ref- 
erence 3 may be used. The 1/7-power law for the velocity 
distribution in the turbulent boundary layer seems to be 
as good as any simple assumption but, of course, cannot be 
relied upon in the presence of a strongly adverse pressure 
gradient or condition approaching those of turbulent sep- 
aration. The equation is as follows: 
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where the subscript c denotes the value of the quantity 
at the transition point. 



CONCLUSION 



The greatest uncertainty about boundary-layer calcu- 
lations such as those outlined is usually associated with 
our lack of knowledge about where transition will occur. A 
criterion suggested, R5 = 8000 to 9500 , may he considered 
only a very rough estimate under favorable conditions. It 
is possible that higher values may be reached in flight. 
Premature transition resulting in lower maximum values of 
Rg may be brought about by air-stream turbulence in wind- 
tunnel work and by vibration and disturbances due to sur- 
face irregularities in flight. To avoid the roughness ef- 
fects, it has been specified that no roughness detectable 
to the finger tips or no waviness that may be felt or that 
might noticeably affect the pressure distribution should 
be allowed, although recent tests have indicated that these 
requirements may be too severe in some cases. 
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